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FOREWORD 

This  report  was  prepared  by  Douglas  L.  Bowers  of  the  Internal 
Aerodynamics  Group,  Aerodynamics  and  Airframe  Branch,  Aeromechanics 
Division  of  the  Air  Force  Flight  Dynamics  Laboratory,  Wriqht-Patterson 
Air  Force  Base,  Ohio.  This  research  was  conducted  under  Work  Unit 
Number  14760221,  "Design  and  Analysis  of  Advanced  Strategic  and  Tactical 
Military  Aircraft  Exhaust  Nozzle  Systems".  This  research  was  conducted 
from  1 January  1976  to  1 September  1976. 
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SECTION  I 
INTRODUCTION 

The  requirement  for  advanced  aircraft  airframe  systems  and  advanced 
propulsion  systems  has  resulted  in  extensive  efforts  in  analytical  and 
experimental  system  development.  In  the  exhaust  nozzle  area,  this  effort 
has  produced  configurations  as  simple  as  the  nozzles  of  the  1945  vintage 
P-80  ''Shooting  Star"  and  the  1949  vintage  F-84  "Thunderjet"  and  as  com- 
plex as  the  blow-in-door  ejector  nozzles  of  the  1964  vintage  F-lll  and 
F-12  and  balanced  beam  nozzle  of  1972  vintage  F-15  and  F-16.  As  exhaust 
nozzle  systems  become  more  complicated,  the  required  wind  tunnel  models 
necessarily  become  more  complex  and  costly.  Limited  by  resources  and 
manufacturing  capability,  wind  tunnel  model  propulsion  systems  are 
simulated  by  high  pressure  cold  air  systems.  This  simulation  has  oro- 
duced  adequate  results  by  neglecting  the  effect  of  a hot  exhaust  gas  com- 
posed of  jet  fuel  combustion  products  on  external  aftbody  nozzle  perform- 
ance. Previous  experimental  work  (References  1,  2)  has  indicated  that 
changes  in  exhaust  gas  temperature  (or  ratio  of  specific  heats)  alter 
the  exhaust  plume  shape,  especially  the  initial  plume  inclination  angle. 
Attempts  to  correct  cold  jet  aftbody  performance  for  hot  exhaust  jet 
effects  using  the  initial  plume  angle  as  a correlating  parameter  have  been 
reasonably  successful  for  a limited  data  base.  This  document  describes 
the  application  of  the  initial  plume  anale  correlation  for  hot  jet  exhaust 
effects  to  an  extended  data  base,  i.e.,  all  available  applicable  hot  exhaust 
jet  data  produced  from  1953  to  the  precent  time.  Also  presented  are  per- 
tinent results  from  other  hot  exhaust  etas  nozzle  test  programs  which  give 
insight  to  the  overall  effect  of  temoorafure  or  ratio  of  snecific  heat 
variations  on  nozzle  boattail  performance.  Described  is  the  general  aftbody 
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flow  problem,  the  methods  of  calculating  the  initial  plume  inclination 
angle,  the  experimental  data,  and  application  of  the  correlation.  Oata 
analysis  was  conducted  from  1 January  76  to  1 September  76. 


I 
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SECTION  II 
DISCUSSION 

1.  BASIC  AFTBODY  NOZZLE  FLOW  PARAMETERS 

The  basic  geometric  parameters  for  nozzle  flow  are  described  in 
Figure  1.  The  nozzle  contour  shown  is  for  a convergent  divergent  nozzle. 

9gy  is  the  trailing  edge  nozzle  external  boattail  angle,  9fJ  is  the  nozzle 
internal  divergence  angle,  and  is  the  initial  plume  inclination  angle 
measured  from  the  nozzle  exit.  The  6j  value  is  defined  for  supersonic 
flow  exhaust  plumes  only.  The  geometric  areas  are  A_,  nozzle  throat  area 
and  Ae>  nozzle  exit  area.  Flow  parameters  are  Pj,  nozzle  throat  total 
pressure,  P , nozzle  internal  static  pressure  at  nozzle  exit,  P-,  static 

“ J 

pressure  at  nozzle  exit  after  internal  expansion,  M , Mach  number  at 

nozzle  exit,  M.,  Mach  number  after  internal  expansion  at  nozzle  exit,  PnT, 

J B i 

boattail  or  base  static  pressure,  P , freestream  static  pressure,  and  M^, 
freestream  Mach  number.  The  ratio  of  specific  heats,  y,  is  subscripted 
"I"  for  the  internal  flow  and  "E"  for  the  freestream  flow. 

For  jets  exiting  into  a subsonic  stream  the  internal  flow  can  assume  four 
conditions;  1)  all  subsonic,  2)  mixed  subsonic/supersonic  flow,  3)  all 
supersonic  with  perfect  expansion  and  4)  all  supersonic  and  underexpanded. 

For  the  all  subsonic  case.  Figure  2a,  (internal  and  external  flow),  the 
static  pressure  at  the  exit  is  matched  internally  and  externally.  In 
Figure  2b,  the  nozzle  total  pressure  is  sufficient  to  have  choked  flow  cl 
the  nozzle  throat  and  M>1  flow  in  a section  of  the  diverging  section.  The 
static  pressure  at  the  exit  cf  the  nozzle  is  matched  internal Iv  and  ex- 
ternally by  a combination  of  the  normal  shock  and  the  subsonic  expansion 
in  the  remaining  divergent  section.  In  condition  2c.  the  pressure  is 
matched  with  a normal  shock  at  the  exit.  In  these  three  cases  2a,  2b,  and 

2c,  the  nozzle  exit  flow  is  similar  and  an  initial  plume  inclination  angle 

3 
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6-  is  not  defined.  The  case  where  the  6.  is  defined  is  shown  in 

J J 

Figure  2d.  The  static  pressure  in  the  nozzle  at  the  exit  after  the 
internal  expansion  is  still  greater  than  the  freestream  static 
pressure.  in  this  case,  the  flow  expands  tc  create  the  plume 
boundary  originating  at  the  nozzle  exit  . When  a distinguishable 
nozzle  plume  exists,  a ij  is  defined.  For  operational  aircraft, 
the  total  oressure  at  the  throat  is  generally  large  enough  to  result 
in  underexpanded  flow  at  the  exit  and  a subseguent  initial  plume  in- 
clination angle.  The  important  parameters  are  Py  (total  pressure  at 
the  throat),  A^/ (nozzle  area  ratio)  and  P^p  (usually  the  pressure 
in  the  freestream,  P , or  PRT,  a static  Dressure  on  the  external 


surface  near  the  exit). 

For  jets  exhausting  into  a supersonic  freestream,  the  flow 
mechanisms  are  shown  in  Figure  3.  When  the  nozzle  throat  total 
pressure  is  not  sufficient  to  provide  M=1  flow  at  the  throat.  Figure 
3a,  the  internal  flow  is  subsonic  throuohout  and  static  pressure  is 
controlled  bv  the  external  reference  pressure  and  matched  at  the 
nozzle  exit.  For  the  mixed  internal  flow  case.  Figure  3b,  the 
internal  normal  shock  location  and  the  amount  of  subsonic  expansion 
•interact  with  the  external  supersonic  flow  after  an  oblique  shock  or 
expansion  system  to  match  static  pressure  and  direction  at  the  nozzle 
exit.  To  solve  this  case  for  the  correct  combination  of  normal  shock 
location  and  external  flow  oblique  shock  or  expansion  system  becomes 
a complicated  problem  dependent  on  throat  total  oressure,  nozzle  area 
ratio,  and  external  reference  pressure.  In  example  3c,  the  flow  is 
supersonic  throughout.  The  matching  of  static  pressure  and  flow 
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Figure  3.  Nozzle  Flow  - Supersonic  Freestream  Mach  Number 
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direction  at  the  nozzle  exit  is  accomplished  by  oblique  shock  or  expan- 
sion wave  systems  and  is  dependent  on  the  external  reference  pressure, 
external  nozzle  boattail  angle,  e^,  internal  nozzle  divergence  angle, 
0^,  and  nozzle  internal  area  ratio.  The  final  matched  flow  direction 
along  the  exhaust  plume  boundary  is  the  initial  plume  inclination  angle, 
6 j . Only  in  case  3c  is  an  exhaust  plume  and  an  initial  plume  inclination 
angle  defined.  Any  effects  of  the  on  the  external  nozzle  pressures 
must  be  transmitted  forward  by  the  nozzle  boundary  layer  . Cases  2d  and 


3c,  where  6,  is  defined,  will  be  emphasized  in  this  document. 

J 

2.  EFFECT  OF  TEMPERATURE  VARIATION  ON  EXHAUST  PLUME  6 . 


Experimental  and  analytical  work  done  in  the  mid-1950‘s  established 
the  effect  of  y,  ratio  of  specific  heats,  on  the  exhaust  plume  <5., . This 
work,  References  3,  4,  and  5,  combined  wind  tunnel  Schelieren  photographs 
and  the  method  of  characteristics  computer  code.  Figure  4,  from  Reference 
4,  shows  the  initial  plume  inclination  angle  increasing  with  decreasing 
y for  a Mj=2.0  and  at  any  given  static  pressure  ratio.  This  variation  of 
<5j  with  y is  a result  of  sensitivity  of  the  Prandtl -Meyer  turning  angle, 
v,  and  static  to  total  pressure  ratio,  (P/Py),  with  y.  In  NACA  1135, 
"Equations,  Tables,  and  Charts  For  Compressible  Flow",  pg  14,  Equations 
171a  and  173c  are  the  following  relations: 


tan 


-1 


'y+1 


(m‘ 


•1))1/2-(90-u 


p rtl  1/2  -1  p 1/2 

(fi— ) Y - J_0  + cos (2  y ■ 1 j (v+  tan  (M  -1)  ))) 

T Y+1  ly+H 
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Figure  4.  (Reference  4)  Initial  Plume  Angle  vs.  Static  Pressure  Ratio 


AFFDL-TR-76-1 38 


These  relations  are  presented  graphically  in  Figures  5 and  6. 

In  hot  jet  exhaust  testing,  the  temperature  is  often  varied  by 
changinq  the  fuel/air  ratio  to  an  ethylene  burner.  This  produces  a 
Y change  as  well  as  a temperature  change  in  the  jet  exhaust.  Temperature 
effects,  therefore,  can  be  approached  as  y effects  and  both  will  be  dis- 
cussed as  hot  exhaust  gas  effects  from  this  point  on. 

3.  EFFECT  OF  EXHAUST  PLUME  6.  ON  AFTBODV  NOZZLE  PERFORMANCE 

In  an  inviscid  flow  field,  the  exhaust  nozzle  plume  is  modelled  as 
a solid  body  interacting  with  the  aftbody  nozzle  flow  field.  The  geo- 
metric character  of  the  plume  at  the  nozzle  exit,  i.e.  6^,  determines  the 
amount  of  turning  the  aftbody  flow  must  negotiate  at  the  nozzle  exit, 
determines  the  aftbody  nozzle  pressure  distribution,  and  hence  the  ex- 
ternal nozzle  performance.  References  6 and  7 indicate  the  boattail 
drag  coefficient  variations  with  the  chanqing  exhaust  plume  6 ^ ' s . 

Figure  7 shows  the  effect  of  nozzle  plume  simulators  on  aftbody  drag 
for  a typical  nozzle.  As  the  solid  plume  simulator  is  chanaed  from  a 
solid  cylindrical  sting  with  a base  area  to  a cylindrical  sting  with 
no  base  area  to  an  underexpanded  plume  shape,  the  6^  increases  and  the 
aftend  drag  coefficient  subsequently  decreases.  Each  of  these  three 
solid  plume  simulators  correspond  to  a nozzle  pressure  ratio  and  a 5; 
value.  The  aftbody  drag  variation  with  <5j  indicates  the  importance  of 
the  exhaust  plume  initial  inclination  angle  on  aftbody  performance. 

4.  CORRECTION  OF  COLD  EXHAUST  GAS  DATA  FOR  HOT  EXHAUST  GAS  EFFECTS 

Typical  advanced  turbofan  and  turbojet  engines  have  exhaust  gas  total 

temperatures  near  3200°R  in  afterburning  power  and  lbOOuR  ir.  dry  cruise 
power.  Past  attempts  to  establish  6j  as  a correlating  parameter  for 
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hot  exhaust  gas  effects  are  reported  in  References  1 and  2.  C,  E. 
Robinson,  ARO,  Inc.  (Reference  1),  investigated  hot  exhaust  gas  effects 
for  temperatures  from  530°R  (cold)  to  approximately  3200°R.  The  effect 
of  the  hot  exhaust  on  total  aftbody  drag  in  this  study  is  shown  in 
Figure  8 for  = 0.9.  As  the  exhaust  gas  temperatures  are  increased  from 
530°R  to  2390°R,  the  aftbody  drag  coefficient  decreases.  Similarly. 

W.  Compton,  NASA  Langley,  Figure  9 (Reference  2),  shows  the  aftbcdy  drag 
coefficient  for  different  y's  and  temperatures  collapsing  together  to  an 
extent  when  correlated  using  5..  The  top  plot  is  versus  Py/P^  or 
nozzle  pressure  ratio  (NPR).  The  middle  plot  is  versus  static  pres- 
sure ratio  and  the  lower  plot  is  Cq  versus  6j . If  the  <5j  correlation 
successfully  accounts  for  all  hot  exhaust  gas  effects,  the  lower  plot  in 
Figure  9 would  consist  of  one  line,  i.e.,  all  data  would  collapse  on  all 
other  data.  It  should  be  noted  that  the  6.  correlation  does  not  entirely 
compensate  for  hot  exhaust  gas  effects  but  does  tend  to  compress  the  data. 


M_  = 0.9 


Figure  8.  (Reference  1)  Drag  Coefficient  Sensitivity  with  Temperature 


14 


AFFDL-TR-76-1 38 


T ,°R  M =0.90 

oo 

540 

1162.3 
1458. 

1323.4 


OPR 


Y 

1 .40 

1.30 

1.23 

1.26 


Figure  9.  (Reference  2)  Aftbody  Drag  Sensitivity 
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To  use  6 . as  a correlating  parameter,  the  following  procedure  is  used. 

J 

(See  Figure  10). 

a.  The  baseline  data  for  a cold  jet,  i.e.  y=1.4  for  air,  is  in  a form 
Cg  vs  NPR  (experimental  data). 

b.  Using  the  techniques  discussed  in  the  next  section,  a plot  of 
6 . vs  NPR  is  made  as  a function  of  y. 

c.  The  vs  NPR  plot  is  entered  at  the  NPR  for  the  data  to  be 
corrected. 

d.  For  a correction  to  y=l . 3 i.e.  a hot  jet,  continue  up  the  plot 
until  the  y=1.3  curve  is  intersected.  This  intersection  gives  the  6.  to 

J 

simulate  using  y=1.4  flow. 

e.  For  this  6 ^ value,  move  laterally  on  the  plot  until  the  y=1.4 
curve  is  intersected.  This  intersection  gives  the  NPR  using  y=1.4 
exhaust  gas  that  simulates  a y=1.3  exhaust  gas,  i.e.  simulates  the  same  6^. 

f.  The  Cp  value  at  this  corrected  NPR  value  is  the  original  drag 
value  corrected  for  hot  exhaust  effects. 

5.  CALCULATION  OF  AND  CONSTRUCTION  OF  THE  6-  VS  NPR  PLOT 

A computational  method  for  will  be  developed  for  the  examples 
shown  in  Figures  2d  and  3c,  i.e.  subsonic  and  supersonic  external  flow 
with  an  underexpanded  plume.  The  technique  described  is  programmed  on 
the  Hewlett-Packard  9830  computer  located  in  the  Aerodynamics  and 
Airframe  Branch,  Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson 
Air  Force  Base,  Ohio.  NACA  1135  equation  numbers  are  listed  in  parenthesis. 

The  calculation  of  the  initial  plume  angle  for  subsonic  external 
flow  is  as  follows: 
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The  known  values  for  this  calculation  must  be  PT ’ Ae’  AT’  PBT’  V 
Yj , and  The  calculation  procedure  is  as  follows: 

1.  Compute  the  area  ratio,  A ^/A-j. 

2.  Using  the  NACA  1135  compressible  flow  tables,  assuming  M=1 
at  the  nozzle  throat,  Ae/A*  indicates  a Mach  number  at  the  exit,  Mg. 

3.  The  exit  static  pressure  is  a function  of  Py,  M , and  -,-j  • 

P =P  (h^Im  )Y'1  (40) 

e T'  2 e' 

4.  The  Prandtl -Meyer  function  of  the  exit  flow  ve,  is  a function 

of  Yt  and  M . 

I e 

V+1  1/2  -I  o 1/2  , ? 1/2 

ve  =(*Tj-)  tan  ((^-)(Mg-l))  -tan  (MM)  (171c) 

5.  The  static  pressure  to  be  matched  along  the  initial  plume  boundary 


is  Pgy.  The  jet  Mach  number' is  a function  of  Pgy  and  Py,  therefore: 
Pr.. 

(44) 


Mi  =(((~-) 
J PT 


1/2 


-ll~) 
f ' 


6.  The  Prandtl-Meyer  function  in  the  jet,  i.e.  after  the  expansion, 
is  determined  as  in  step  4 with  v,  as  a function  of  Mi. 

J J 

vj  = FfMj.y)  (171c) 

7.  For  the  Prandtl -Meyer  expansion, 

Vve^V  vj 

where  0 and  0.  are  flow  directions  with  resnect  to  the  nozzle  sidewall, 
e j 

0p,  therefore,  is  zero.  0j  is  of  negative  sign  if  the  flow  direction  is 
away  from  the  nozzle  center! -'ne.  A new  variable  8. a is  ot  opposite  sign 
convent  ion , i .e.  positive  if  the  flow  direction  is  away  from  trie  nozzle 
center 1 in?  , Figure  ' 1 . 
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Therefore, 


0.  = v„  - v. 

J e J 

0 .+  = v . - v 
J J e 

8.  To  convert  to  , the  flow  direction  must  be  referenced  to  the 
horizontal.  The  correction  to  Qj+  is  the  nozzle  divergence  half  angle. 

6j  = °j+  + °N 

To  construct  the  6'j  vs.  N PR  plot  using  this  technique,  the  nozzle 

throat  total  pressure  is  varied,  corresponding  to  the  experimental 

nozzle  pressure  ratio  values  and  the  subsequent  6-  is  read  from  the 

J 

compute*"  solution. 

The  supersonic  external  flow  6j  calculation  process  involves 
matching  the  flow  direction  and  pressure  at  the  exit  for  the  'eternal 
and  external  nozzle  flow  fields.  The  internal  flow  at  the  exit  can  be 
adjusted  by  either  an  oblique  shock  or  a Prandt-rleyer  expansion.  Since 
6.-  is  defined  only  for  angles  above  the  horizontal  plane,  the  external 

O 

flow  changes  are  related  to  oblique  shock  solutions  only.  To  achieve 

the  proper  match  at  the  exit,  a plot  of  P.  versus  6.  is  constructed. 

J J 

The  jet  pressure,  P • , after  the  Praadtl -Meyer  expansion,  is  calculated 

for  both  internal  and  external  flows  for  a series  of  6^  values.  The 

point  where  the  internal  and  external  flow  curves  cross  is  the  point 

where  P.  and  6-  are  matched  for  the  two  flows.  This  is  the  correct  6. 

J J J 

value  for  a given  set  of  flow  conditions. 
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For  the  internal  flow  solution,  the  calculation  proceeds  as 
follows : 

a.  Compute  the  area  ratio,  Ae/Ay. 

b.  Using  compressible  flow  tables,  assuming  M=1  the  nozzle 
throat,  Ag/A*  provides  an  exit  Mach  number,  Mg. 

c.  The  exit  static  pressure  is  a function  of  P^,  Mg  and 

P =p  (i+lli  M 2):^T  (40) 

6 1 2 e 

d.  The  oblique  shock  solution  is  started  for  a wedge  angle 
equivalent  tc  the  9^  value.  This  gives  a 5.-  of  0.  (Figure  12). 

e.  From  NACA  1135,  Equation  151,  the  static  pressure  ratio 
across  the  oblique  shock  is  approximated  by: 


P,  yM26 

Pe  /„2 


(MM) 


' T72+Y' 1 


2(y^1)M4-^(M2-1),2 
4(M2-!)2 


( 2il)  2m8_ Zil^3v2MG 
+ „ HJ7'  32  n 24 

(m2-i  ) 


+|(y+l)M4-M2+-2-)63+---' 


(151) 


where  M=Mg,  y=>j,  and  6 is  the  angle  through  v.'hich  the  flow  is  turned 
( i n rad  i an  s' . 

f.  The  jet  pressure  is  then 

P^Pct  — ] (151) 
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q.  The  flow  turning  angle,  <5,  is  decreased  by  1 degree  and  the 
calculation  in  steps  5 and  6 repeated.  This  procedure  continues  up  to  and 
including  6=0.  At  that  point,  6.=9fJ.  To  increase  6.,  the  flow  must 
expand  around  the  nozzle  exit. 

h.  The  Prandtl -Meyer  expansion  procedure  is  started  with  an 
expansion  angle  of  0 degrees.  This  point  serves  as  a check  for  the 
last  oblique  shock  point. 

i.  The  Prandtl-Meyer  function  at  the  exit,  v , is  a function  of 
Yj  and  Mg. 

V+1  1/2  , , 2 1/2  , 2 1/2 

tan  '!))  ‘tan  (Me  -1)  (171c) 

j.  From  Method  of  Characteristic  Theory. 


9 +\j  =0  .+v  • 
e e j vj 

where  0 and  0 are  flow  directions  with  respect  to  the  nozzle  sidewall 
® J 

0 , therefore,  is  zero.  0.  corresponds  to  6,  or  flow  turning  angle. 

® J 

Therefore 


TV6 

where  6+  is  positive  if  the  flow  direction  is  away  from  the  nozzle 
centerl ine. 

k.  The  initial  plume  angle,  6 . , is  defined  as: 

W6+ 

l.  The  jet  Mach  number  is  a function  of  v-  and  y,.  Use 
compressible  flow  tables  to  get  M ■ . M.-F(vj,  Yj) 

The  jet  static  pressure  is  a function  of  P^,  M;,  and  y,. 


m. 


-V 
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n.  The  turning  angle,  6,  is  incremented  by  2 degrees  and 

steps  10  through  13  are  repeated.  The  calculation  is  arbitrarily 
stopped  when  the  flow  has  expanded  through  the  0^+10  degrees. 

The  external  flow  solution  is  an  oblique  shock  computation 
only.  For  an  initial  plume  inclination  angle  of  0,  the  flow  must  turn 
back  through  an  angle  equivalent  to  the  external  boattail  angle,  0^^. 

The  calculation  is  started  at  that  point. 

a.  The  initial  turning  angle,  6,  must  be  equivalent  to  6gy,  Figure  13. 

6j=0;  6=0bt 

b.  As  in  step  5 of  the  internal  flow  oblique  shock  computation. 
Equation  151  in  MACA  1135  approximates  the  pressure  jump  across  the 
oblique  shock. 

Pj=PT  [-(151)-] 

oo 

c.  The  turning  angle,  6,  is  incremented  by  one  degree  and  the 
calculation  is  continued  at  step  2.  The  calculation  is  arbitrarily 
completed  wi  th  5-  (0^+10°),  since  a pressure/d  ^flection  match  has  usually 
occurred. 

From  these  calculations  for  external  and  internal  supersonic 
flow,  the  P.  versus  <S . plot  is  drawn.  The  intersection  of  the  two  curves 

J J 

gives  the  flow  direction  and  pressure  at  the  nozzle  exit  for  the  given 

flow  conditions,  NPR,  PBT  ’ etc.  A sample  calculation,  which  has  been 

verified  by  hand  calculations,  is  presented  in  Figure  14.  The  P..  versus 

%> 

plot  in  Figure  15  is  a result  of  this  sample  calculation. 
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Figure  14.  Sample  Supersonic  Freestream  Calculation 
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Figure  15.  Pressure/Deflection  Match  Point 
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6.  SURVEY  OF  AVAILABLE  EXPERIMENTAL  DATA  APPLICABLE  TO  HOT 

EXHAUST  GAS  EFFECTS 

This  section  contains  experimental  data  which  is  directly  applicable 
to  an  initial  plume  correlation.  Also  included  is  experimental  data 
which  assesses  the  effect  of  temperature  ana  y on  other  aspects  of  aftbody 
nozzle  testing  techniques. 

With  aftbody  drag  coefficient  plotted  versus  , a successful  initial 
plume  angle  correlation  compresses  all  data  from  different  temperature 
exhaust  gases  onto  one  line.  The  best  available  data  is  found  in  Ref- 
erences 1,  2,  8,  and  9.  C.  E.  Robinson,  et.  al . , Reference  1,  performed 
a combined  experimental  and  analytical  research  program  on  an  isolated 
nozzle  aftbody  wnich  was  tested  at  Mach  numbers  from  0.6  to  1.5  in  the 
Arnold  Engineering  Development  Center's  (AEDC)  transonic  Propulsion  Wind 
Tunnel  (16T).  The  experimental  data  included  static  pressure  orifices, 
force  balance  outputs,  and  Schlieren  photographs.  Both  air-cooled  and 
water-cooled  configurations  were  tested  but  die  to  the  influence  of  the 
secondary  cooling  air  on  the  aftbody  pressure  distribution,  only  the 
water-cooled  configuration  will  be  discussed. 

An  ethylene  burner,  with  combustion  products  similar  to  those  of 
JP-4,  provided  the  model  with  exhaust  gas  temperatures  up  to  approxi- 
mately 2500°R.  Model  details  are  shown  in  Figures  16  and  17. 

The  effect  of  increased  exhaust  gas  temperature  is  shown  in  Figure  18. 

As  the  exhaust  gas  temperature  increases,  the  aftbody  nozzle  drag  decreases. 
Figure  19  is  a plot  of  pressure  coefficient  distribution  versus  axial 
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Flameholder  Flow  Mixer 


Figure  17  (Reference  1)  Model  Detail 
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distance  as  a function  of  exhaust  gas  temperature.  A temperature  increase 
raises  the  pressure  coefficient  values  on  most  of  the  afterbody  and  nozzle 
contours.  If  this  change  in  pressure  coefficient  distribution  can  be 
attributed  to  a change  in  the  initial  plume  angle  the  data  should  collapse 
onto  one  line  when  aftbody  drag  coefficient  is  plotted  versus  6j.  This 
correlation  is  shown  in  Figure  20.  The  correlation  does  tend  to  draw  the 
data  together,  but  still  does  not  compensate  for  the  2000CR  case  and  tends 
to  over-correct  the  1420°R  data  at  hioher  nozzle  pressure  ratios. 

At  a freestream  Mach  number  of  1.5,  the  effect  of  exhaust  gas  tempera- 
ture is  small  (Figure  21).  When  the  drag  coefficient  data  is  plotted 
versus  the  initial  plume  angle,  the  effect  of  the  hot  exhaust  gas  is  over- 
corrected for  the  middle  range  temperatures  (Figure  22). 

Additional  hot  jet  data  from  the  Arnold  Engineering  Development 
Center's  16T  facility  is  described  in  Reference  g.  The  test  apparatus 
(Figure  16)  is  identical  to  that  used  by  Robinson  (Reference  1)  with  the 
exception  of  the  nozzle  boattail  contours.  These  aftbody  boattails, 
designated  the  15°  and  25°  AGARD  boattails,  .ere  tested  b.v  other  MATO 
countries  under  the  AGARD  program  "Improved  Nozzle  Testinq  Techniques  in 
Transonic  Flow".  Efforts  of  other  NATO  countr'es  will  be  discussed  in 
additional  paragraphs  of  this  report. 

Using  the  initial  plume  angles  calculated  by  AEDC,  the  measured 
aftbody  drag  coefficients  are  plotted  versus  nozzle  pressure  ratio  and 
initial  plume  angle  in  Figures  23  through  27.  The  figures  include 
aftbody  drag  coefficients  versus  nozzle  pressure  ratio  and  initial  plume 
inclination  anale  for  the  15°  AGARD  boattail  at  Mach  numbers  of  0.8,  0.9, 
and  1.1.  In  all  cases,  an  increase  of  exhaust  nas  temperature  resulted 
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Figure  21.  Drag  Coefficient  vs.  NPR,  Reference  1 Data,  M =1.5 
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in  a decrease  in  aftbody  drag.  The  exhaust  gas  temperatures  listed  are 
nominal  values.  When  the  drag  coefficient  is  plotted  versus  the  initial 
plume  angle,  the  plots  of  different  temperatures  tend  to  collapse  on  one 
another,  especially  at  the  higher  values  of  6j  , corresponding  to  the  higher 
nozzle  pressure  ratios.  At  low  values  of  dL  , the  correlation  does  not 
appear  to  be  effective.  It  should  be  noted  that  the  operating  nozzle  pres- 
sure ratios  of  current  turbofan  and  turbojet  engines  in  the  cruise  con- 
ditions is  3 to  6.  At  this  area  of  interest,  therefore,  the  6^  correlation 
is  not  sufficient  to  handle  hot  exhaust  effects. 

Figures  26  and  27  present  drag  coefficient  data  for  the  25°  AGARD  boat- 
tail.  This  configuration  has  a region  of  highly  separated  flow  from  the 
boattail  shoulder  aft  to  the  exit.  The  aftbody  drag  coefficient,  therefore, 
is  relatively  insensitive  to  nozzle  pressure  ratio  changes.  When  the  drag 
coefficient  is  plotted  versus  initial  plume  angle,  the  data  spread  due  to 
hot  exhaust  gas  effects  is  unchanged.  This  data  plus  the  low  nozzle  pres- 
sure ratio  data  from 'the  15°  AGARD  boattail  indicates  that  areas  where  the 
drag  coefficient  is  marginally  influenced  by  nozzle  pressure  ratio  are  areas 
where  the  initial  plume  angle  correlation  is  ineffective. 

The  aftbody  drag  trends  and  the  applicability  of  the  initial  plume  angle 
correlation  to  a hot  exhaust  jet  are  substantiated  by  Reference  2.  VI.  R. 
Compton,  using  an  isolated  nacelle  with  cold  air  or  hydrogen  peroxide  ex- 
haust jets,  investigated  the  magnitude  of  the  effect  of  various  jet  parameters 
on  external  nozzle  boattail  performance. 
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Figure  27.  Drag  Coefficient  vs.  NPR  and  Initial  Plume  Angle 

Correlation,  Referem  e 9 Rata.  PS  Coattail,  M 0.9 
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The  test  apparatus  and  nozzle  configurations  are  shown  in  Figure  28. 
Two  nozzles  with  external  boattail  angles  of  10°  and  20°  were  tested  with 
design  internal  area  ratios  to  provide  Mp=l  and  Mg=2.  Primary  instrumen- 
tation was  aftbody  external  pressures,  surface  temperatures,  nozzle  exit 
pressures  and  temperatures,  and  external  boattail  boundary  layer  profiles 


Figure  28.  (Reference  2)  Nozzle  Details 
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As  the  exhaust  gas  temperature  increased,  the  aftbody  nozzle  pressure 
distribution  was  altered  to  generally  give  higher  pressures.  Figure  29 
shows  this  trend  for  a freestream  Mach  number  of  0.9.  This  altered 
pressure  distribution  produces  the  reduced  nozzle  drag  coefficient.  One 
of  the  results  of  this  investigation  is  that  the  greatest  drag  differences 
for  varying  exhaust  gas  temperature  occurred  for  high  transonic  Mach 
numbers  (M»0.95)  and  large  external  boattail  angles  (6gy=20°).  A 
correlation  of  the  drag  coefficient  and  the  initial  plume  angle  is  shown 
in  Figures  30  and  31  for  freestream  Mach  numbers  of  0.9  and  1.2, 
respectively.  The  correlation  tends  to  collapse  the  drag  coefficient 
data  together,  especially  at  large  values  of  nozzle  pressure  ratio. 

This  trend  is  consistent  with  previous  work.  Compton  attributed  the 
remaining  spread  in  drag  coefficient  data  to  jet  entrainment. 

Tests  conducted  by  the  NACA  on  simple  convergent  nozzles  demonstrated 
the  effects  of  a hot  exhaust  jet  on  afterbody  nozzle  drag.  The  transonic 
tests  were  conducted  in  the  Langley  8-foot  transonic  wind  tunnel  with 
nacelle  contours  of  varying  boattail  angles  (Figure  32).  Exhaust  jet 
temperatures  of  540°R,  1260°R,  and  1660°R  were  provided  by  an  ethylene 
burner.  Primary  instrumentation  was  external  aftbody  static  pressures. 

Using  the  computer  coding  for  subsonic  external  flow,  the  initial 
plume  angle  correlation  was  computed  at  Mach  numbers  of  0.8  and  0.9.  Plots 
of  total  aftbody  drag  coefficient  versus  nozzle  pressure  ratio  (NPR)  and 
initial  plume  inclination  angle  (iS^)  are  presented  in  Figures  33  and  34 
for  Mach  numbers  of  0.8  and  0.9,  respectively.  For  aftbody  boattail  con- 
figurations with  relatively  low  boattail  angles,  i.e.  16  and  24  degrees. 
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Figure  33.  (Reference  8)  Drag  Coefficient  vs.  NPR  and  Initial  Plume  Angle 
Correlation,  M =0.8 
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the  initial  pi  time  angle  correlation  tends  to  collapse  the  data  together 
correcting  for  most  if  not  all  effects  of  the  hat  exhaust  gas.  For  boat- 
tails  with  large  boattail  angles,  i.e.  0BT=45O,  the  initial  plume  angle 
correlation  tends  to  leave  the  data  spread  unchanged  or  increase  the  data 
spread. 

In  summary,  the  initial  plume  angle  correlation  apoears  to  be  effective 
for  high  nozzle  pressure  ratios  and  moderate  temperature  increases.  For 
high  temperature  exhaust  flows,  separated  flow  boattails,  or  attached 
flow  boattails  at  low  nozzle  pressure  ratios,  the  initial  plume  angle 
correlation  does  not  correct  for  hot  exhaust  gas  effects  and  in  some 
cases,  over-corrects  the  drag  coefficient  data.  This  trend  is  supported 
in  part  by  the  four  data  sources  listed  above. 

While  the  amount  of  hot  exhaust  qas  experimental  data  directly  applicable 
to  an  initial  plume  angle  analysis  is  limited,  numerous  other  related  sources 
exist  which  provide  trends  and  information  pertaining  to  the  effect  of  a hot 
exhaust  jet  on  nozzle  boattail  performance. 

Ratio  of  specific  heat  chanqes  as  well  as  temperature  changes  are  pro- 
duced by  a air/fuel  combustion  process.  Py  using  different  gases  in  an 
exhaust  jet,  the  v variation  can  be  produced  without  the  accompanying 
temperature  rise.  Early  efforts  by  N.'.CA  and  NASA,  References  10,  II,  and 
12,  utilized  cold  air  (y=1.4,  T=540°R),  heated  air  (y=1.3,  T=2100"R), 
and  a cold  carbon  dioxide  jet  (y=1.3,  T=540°R).  Lee,  Reference  10, 
investigated  jet  plume  shape  and  cold  jet  simulation  of  hot  exhaust 
jets  at  transonic  Mach  numbers  at  the  NASA  Ames  Research  Center. 

Primary  instrumentation  for  the  strut  supported  cone  cylinder  model 
with  a convergent  nozzle  was  Schlieren  photoqraphs  and  flow  field  pressure 
surveys . 
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Concerning  the  jet  plume  shape,  the  test  results  showed  that  for  the 
hot  air  jet  and  the  carbon  dioxide  jet,  both  with  7=1.29,  the  jet  plume 
shape  was  similar  at  the  nozzle  exit  and  over  most  of  one  exhaust  plume 
wavelength.  The  cold  exhaust  jet,  7=1.4,  could  be  made  to  match  the  hot 
jet  shape  at  the  nozzle  exit  if  the  cold  jet  nozzle  pressure  ratio  was 
increased.  Downstream,  however,  the  cold  jet  at  an  increased  nozzle  pres- 
sure ratio,  did  not  match  the  hot  jet  plume  shape.  The  result  is  that 
simulation  of  the  hot  jet  boundary  for  at  least  one  wavelength  with  a cold 
jet  must  be  done  at  the  same  7 and  nozzle  pressure  ratio.  It  was  also 
found  that  the  external  shock  wave  location  of  the  hot  air  jet  could  be 
matched  by  both  the  carbon  dioxide  jet  and  the  cold  exhaust  jet. 

The  flow  field  static  pressure  surveys  and  the  resulting  forces  and 
moments  due  to  these  pressure  distributions  indicated  that  the  effect  of 
temperature  was  small  but  the  effect  of  the  ratio  of  specific  heats  change 
was  large.  The  hot  air  data  and  the  carbon  dioxide  data  are  similar  (same 
7,  different  gas  temperatures)  while  the  cold  air  jet  data  is  decidedly 
different  (same  gas  temperature  as  carbon  dioxide,  different  >). 

The  most  pertinent  result  from  this  investigation  is  that  the  effect 
of  a hot  exhaust  jet  flow  pattern  is  simulated  by  a cold  jet  if  the  nozzle 
pressure  ratio  and  ratio  of  specific  heats  can  be  matcned. 

References  11  and  12  also  investigated  gases  of  different  y's  and  the 
same  temperature  including  heated  air.  Rase  and  boattail  pressures  were 
measured  on  the  model  shown  in  Figure  35.  The  chanoe  in  the  bast  pressure 
coefficient  for  two  gases  at  the  same  temperature  but  different  >'s  is 
shown  in  Figure  36.  Also  shown  is  the  effect  on  tho  base  pressure  co- 
efficient of  temperature  and  7 variations  using  air  as  a lot . in  both 
cases, i.e.  a smaller  7 far  the  carbon  dioAido  >m  than  the  cold  air  jet 
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and  the  temperature  and  y variations  for  a heated  and  cold  air  jet,  the 
results  were  to  increase  the  base  pressure  coefficient  and  therefore  de- 
crease the  total  aftbody  drag. 

Correcting  the  y = 1.4  base  pressure  coefficients  for  a gas  with  a 
different  y value  is  done  by  adjusting  the  nozzle  pressure  ratio  to  give 
the  same  initial  plume  angle.  These  corrected  base  pressure  coefficients 
are  shown  as  dashed  lines  in  Figure  36.  This  correction  appears  to  be 
satisfactory  down  to  y values  of  1.3.  This  result  is  consistent  with 
earlier  trends  where  the  correction  to  the  same  initial  plume  angle  was 
good  for  moderate  temperature  changes  and  high  nozzle  pressure  ratios. 

A quanitative  analysis  of  the  effect  of  a heated  jet  on  the  exhaust 
plume  spreading  was  conducted  by  Rousso  and  Kcchendorfer  in  Reference  13. 
Using  a total  pressure  and  total  temperature  probe,  the  exhaust  plume 
downstream  of  a convergent  nozzle  was  surveyed.  Figure  37.  The  hot  exhaust 
jet  was  provided  by  burning  propane  with  the  source  plenum  air.  The 
primary  chamber  temperatures  are  540°R  and  l‘'00°R  with  corresponding 
ratio  of  specific  heat  values  of  1.4  and  1.35. 

Figure  38  shows  how  the  Mach  number  distribution  within  the  plume 
varies  for  a cold  and  a hot  exhaust  jet.  M is  the  local  Mach  number  as 
determined  by  the  probe  and  is  the  calcula-ed  Mach  number  assuming  an 
isentropic  expansion  from  the  nozzle  internal  pressure  to  the  freestream 
pressure.  A value  of  M/M-=0.11  is  assumed  to  be  the  outer  plume  boundary. 
The  Mach  number  profiles  for  the  hot  exhaust  jet  are  similar  to  the 
cold  jet  in  shape  but  show  higher  Mach  numbers  near  the  nozzle 
centerline.  The  effect  of  the  hot  jet  on  the  Mach  number  profiles 
appears  to  diminish  at  downstream  positions  of  the  exhaust  plume. 
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Figure  38  (Reference  13)  Mach  Number  Profiles  for  Cold  and  Hot 
Exhaust  Jets 
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An  additional  area  For  concern  in  hot  exhaust  testing  of  nozzle 
afterbody  configurations  is  external  nozzle  surface  heating.  References 

14  and  15,  conducted  for  AGARD  by  German  and  Netherland  research  agencies 
respectively,  shojv  that  uncooled  nozzle  aftbody  models  can  produce 
anomalies  when  investigating  the  effect  of  a hot  exhaust  jet  on  the 
nozzle  drag  coefficient. 

In  Reference  14,  the  model  is  a HFB  320,  a commercial  transport 
nacelle.  Primary  instrumentation  is  external  static  pressure  orifices 
and  surface  thermocouples.  Exhaust  gas  temperatures  range  from  522°R  to 
1422°R.  The  variation  of  nozzle  drag  coefficient  with  Mach  number  and 
temperature  is  shown  in  Figure  39.  The  higher  nozzle  drag  coefficient 
for  the  hot  jet  case  than  the  cold  jet  is  opposite  of  trends  experienced 
in  other  test  facilities.  The  anomaly  appears  to  be  a result  of  body 
heating  by  tne  exhaust  jet.  The  body  surface  temperature,  recorded  by 
the  surface  thermocouples,  is  approximately  the  jet  exhaust  temperature 
of  1422°R.  This  model  boundary  layer  temperature  gradient  creates  an 
outward  flow  from  the  model.  This  modifys  the  boundary  layer  profiles 
to  a more  turbulent  shape  and  effectively  creates  a "olcwn"  boundary  layer. 
The  more  turbulent  boundary  layer  for  the  heated  body  results  in  less  re- 
compression (Figure  40)  on  the  aftbody  and  therefore  more  drag  for  a hot 
jet  than  a cold  jet. 

Reference  15  also  snows  the  effect  of  body  heating  tends  to  reverse 
the  previous  trends  of  nozzle  aftbody  drag  coefficient  with  exhaust  gas 
temperature.  For  this  data,  however,  the  anomaly  does  not  exist  above- 
nozzle  pressure  ratios  of  3,  Figure  41. 

Both  of  those  references  indicate  that  as  a testing  technique,  care 
must  be  taken  to  avoid  coirpromi si nq  the  external  nozzle  pressure  data  by 
nozzle  surface  heating. 
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Figure  40.  (Reference  14)  External  Nozzle  Aftbody  Pressure  Coefficient 
Distribution  Variations  with  Exhaust  Plume  Temperature 


65 


AFFDt-TR-76-1 38 


SECTION  III 
SUMMARY 

The  application  of  the  initial  plume  angle  as  a correlatino  factor 
for  hot  exhaust  gas  effects  for  nozzle/aftbody  pressure  drag  has  been 
investigated  over  an  extensive  data  base.  In  addition,  general  testing 
techniques  and  flow  characteristics  for  hot  exhaust  gas  simulation  have 
been  summarized.  The  summary  of  results  of  this  effort  are  as  follows: 

1.  As  the  exhaust  gas  temperature  increases,  the  aftbody  nozzle 
pressure  drag  decreases  for  cooled  models.  For  uncooled  models  where 
the  model  surface  temperature  increases  with  increasing  exhaust  gas 
temperature,  the  trend  with  exhaust  gas  temperature  is  to  increase  the 
nozzle  aftbody  drag.  This  is  a result  of  temperature  changes  in  the 
boundary  layer  and  its  effect  on  the  flow  pattern. 

2.  The  initial  plume  angle  correlation  successfully  accounts  for 
hot  exhaust  gas  effects  on  nozzle  aftbody  drag  if: 

a.  the  aftbody  boattail  angle  is  lot  to  moderate,  i.e.  up  to 

15°  and 

b.  the  nozzle  pressure  ratio  is  moderate  to  high,  i.e.  NPR>6,  and 

c.  the  temperature  range  is  not  larr.e  between  the  cold  flow 
exhaust  gas  and  the  hottest  exhaust  qas,  i.e.  540"R  and  2000°R. 

3.  The  initial  plume  anale  correlation  fails  to  account  for  hot 
exhaust  gas  effects  on  nozzle  aftbody  drag  if: 

a.  the  aftbody  has  a large  boattail  angle,  ^BT=^*  or  creator, 
where  the  external  flow  is  senarated  and  nozzle  pressure  drag  is  not 
sensitive  to  changes  in  nozzle  pressure  ratio,  or 
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b.  the  nozzle  pressure  ratio  is  below  moderate  values,  i.e. 
NPR<6, (operati ng  nozzle  pressure  ratios  for  existing  turbojet/turbofan 
engines  is  approximately  2.5  to  6)  or 

c.  the  temperature  difference  is  large  between  the  cold  jet  and 
the  highest  exhaust  gas  temperature,  e.g.  if  AT>2000°R. 

4.  For  aftbodies  where  the  exhaust  qas  temperature  effect  on  aftbody 
drag  is  small,  the  initial  plume  angle  correlation  may  over-correct  for 
the  hot  gas  drag  coefficient  data. 

5.  Simulation  of  the  hot  plume  boundary  for  at  least  one  plume  wave 
length  may  be  made  by  a cold  gas  at  the  same  y and  nozzle  pressure  ratio 
as  the  hot  exhaust  gas. 

6.  The  initial  plume  angle  correlation  treats  the  plume  as  a solid 
body.  The  correlation  fails  in  areas  where  jet  entrainment  and  mixing 
are  important  parameters.  The  initial  plume  angle  correlation  when  com- 
bined with  a correlating  factor  for  jet  entrainment  and  mixing  may  offer 
a better  technique  for  correcting  cold  jet  aftbody  pressure  drag  for  hot 
exhaust  oas  effects. 
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